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ABSTRACT 

We have extended the analysis of lJaffe et al.1 (2005a. 2006) to a complete Markov chain Monte 
Carlo (MCMC) parameter space study of the Bianchi type Vllh models including a dark 
energy density, using Wilkinson Microwave Anisotropy Probe (WMAP) cosmic microwave 
background (CMB) data from the 1 -year and 3-year releases. Since we perform the analysis in 
a Bayesian framework our entire inference is contained in the multidimensional posterior dis- 
tribution from which we can extract marginalised parameter constraints and the comparative 
Bayesian evidence. Treating the left-handed Bianchi CMB anisotropy as a template centred 
upon the 'cold-spot' in the southern hemisphere, the parameter estimates derived for the total 
energy density, 'tightness' and vorticity from 3-year data are found to be: fl to t — 0.43 ± 0.04, 
h = 0.32±g;?|, w = 9.7±{\l x lO" 10 with orientation 7 = 337°±23)- This template is pre- 
ferred by a factor of roughly unity in log-evidence over a concordance cosmology alone. A 
Bianchi type template is supported by the data only if its position on the sky is heavily re- 
stricted. All other Bianchi Vllh templates including all right handed models, are disfavoured. 
The low total energy density of the preferred template, implies a geometry that is incompati- 
ble with cosmologies inferred from recent CMB observations. Jaffe et ail (|2005b) found that 
extending the Bianchi model to include a term in Qa creates a degeneracy in the fi TO — J7a 
plane. We explore this region fully by MCMC and find that the degenerate likelihood con- 
tours do not intersect areas of parameter space that 1 or 3 year WMAP data would prefer at 
any significance above 2a. Thus we can confirm the conclusion that a physical Bianchi Vllh 
model is not responsible for this signature, which we have treated in our analysis as merely a 
template. 

Key words: cosmological parameters - cosmology:observations - cosmology:theory - cos- 
mic microwave background 



1 INTRODUCTION 

It has recently been proposed that some of the current anoma- 
lous results seen in observations of the CMB sky by WMAP, 
namely non-Gaussian structure ce ntered on the so called 'cold 
spot' in the southern hemisphere iVielva et al.l 1 20051 ; ICruz et all 
2005), low quadrupole Efstathiou 2004) and multipole alignments 
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20031) . could be removed by allowing for 



large scale vorticity and sh ear. Components su ch as these arise in 
Bianchi type Vllh models ( Bar row et al Jl 19851) due to large scale 
rotation of the universe distorting the CMB since the time of last 
scattering. This assumption has the unfortunate side-effect of vio- 
lating universal isotropy and hence the cosmological principle, so 
any such result should be examined very carefully. 

Such a finding was indeed made by I Jaffe et all (2005a) [here- 
after referred to as Jaffe] where a significant correlation was found 
between a class of Bianchi Vllh models and WMAP observations. 



However t he resultant best-fit B ianchi component was derived di- 
rectly from Barr ow et al.l il985l) for a universe with a large negative 
curvature with total energy density Sltot = 0.5 and including no 
dark energy component. Such a cosmology cannot be reconciled 
with cosmologies inferre d from either curren t CMB or other astro- 
nomical observations. In I Jaffe et alj d2005bl) the authors extended 
their study to include a cosmological constant and searched for a 
morphologically identical template to their first analysis in the in- 
creased parameter space. We perform a new study to explore fully 
the parameter space of both Bianchi and cosmological parameters 
in order to select the best fitting template. We carry out this in- 
vestigation using Markov chain Monte Carlo (MCMC) sampling 
in a fully Bayesian manner. Hence we also find the comparative 
Bayesian evidence allowing an assessment of whether the inclusion 
of such a template is supported by the data. 
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2 MCMC ANALYSIS AND MODEL SELECTION 
FRAMEWORK 

Several recent studies have been performed using a Bayesian ap- 
proach to cosmological parameter estimation and model selection 
djaffdll996l iDrell et alJl200Ct IJohn & Narlikajliijoi, Isiosar et al ' 
2003L ISaini etalJ |2004 [Marshall et al.l l2003l iNiarchou et al 
2004 iBasset et al J 1200-4 iMukheriee et alj 12001 iTrottal l2o6i 



Beltran et alj d2005h 7 lBridges et al.l I2005L IParkinson et alj |2006h . 
Here we are interested in determining the best fitting Bianchi Vllh 
sky from WMAP data using a selection of all sky maps avail- 
able; the 1-year WM AP Internal Linear Combination (WILC) map 
jBennett et alj|2003[). the Lagran ge Internal Linear Combination 
(LILC) map dEriksen et al. 2004 ), the Tegmar k, de Oliveira-Costa 
& Hamilton (TOH) map (Tegma rk et alj|2003l). the Maximum En- 
tropy Method (ME M) map dStolvarov et alj l2005t) and the new 
3-year WILC map dHinshaw et alj l2006h . We have incorporated 
Bianchi s ky simulations into an adapted version of the COSMOMC 
package dLewis & Bridlell2002h which we show via simulated sky 
maps is capable of extracting both the underlying cosmology and 
the Bianchi signal, while also determining when the addition of 
such a component is favoured by the data. 

The Bayesian framework allows the examination of a set of 
data D under the assumption of a model M defined by a set of 
parameters © with B ayes' Theorem: 



P(0|D,M) 



P(D|Q,M)P(Q|M) 
P(DIM) ! 



(1) 



where P(0|D, M) is the posterior distribution, P(D ©, M) the 
likelihood, P(©|M) the prior and P(DjM) the Bayesian evi- 
dence. We explore the parameter space by MCMC sampling to 
maximise the likelihood utilising a reas onable convergence crit eria 
from the Gelman & Rubin R statistic dGelman & Rubinlll992h so 
that once the Markov chain has reached its stationary distribution, 
samples from it reliably reflect the actual posterior distribution of 
the model. The final result will be parameter estimates marginalised 
from the multi-dimensional posterior and a value of the evidence 
for the model being studied. 

2.1 Likelihood 

A CMB experiment records the temperature anisotropy on the sky 
T(p). We wish to find a cosmological model that can optimally de- 
scribe these fluctuations by maximising the probability of observ- 
ing the data T(p) given a model described by a set of parameters 
©c- Then for a Gaussian model, from Bayes' Theorem we can ex- 
press this probabilty, the likelihood, as: 



P(T(p)|0c) 



5 -|(T T C~ 1 T) 



(2) 



where C is the covariance matrix containing the predictions of 
the model ©c and \C | its determinant. Expanding the temperature 
sky in spherical harmonmics T(p) = 5_a m a imYi m , in a uni- 
verse that is globally isotropic means Ci m ji m i = (ai m a*, m ,) = 
CiSuiS mrn i where the variance of the harmonic coefficients ai m is 

C* — WpL Em l a '« 



| 2 . Thus Eqn.[2]becomes0: 

e (-\a lm \ 2 /2C t ) 



P({ai m }\®c) oc H ■ 



im 



(3) 



1 Where we have used the notation that C; and af m refer to observed quan- 
tities while Ci and a; m refer to their theoretical counterparts. 



If one assumes universal isotropy the likelihood is then independent 
of m which can be summed over to give (for a full sky survey) the 
form used in most conventional CMB analyses: 



In Pi C ';!©,■) = (2/ + I ; ( gi+ln 



c, 



(4) 



where Ci is the observed quantity. However Bianchi Vllh models 
are anisotropic so one cannot simply compare observed Ci 's with 
theoretical C'i's, instead the full set of a; m 's must be used. 

In this analysis we make the assumption that there is both a 
background cosmological contribution and Bianchi Vllh contribu- 
tion to the CMB. We take the cosmological component to be that 
of a standard ACDM universe, producing a set of a; m 's, while the 
Bianchi model component (defined by parameters ©b) is given by 
af m . The Bianchi component can then be subtracted from the ob- 
served aim's so that Eqn[3]becomes: 



P({a !m }|© B ,e c )c< n^F e 



/ i - B ,2 N 

m — 1 

For numerical convenience we take the natural logarithm to yield 
the final likelihood function used in this analysis, hereafter referred 
to as C: 



(aio - af f 
2Ci 



— ^2 \\ aim ~ a fm\ 2 



C, 



(6) 



which would reduce to Eqn.|4]if the Bianchi signal were zero. 



2.2 Evidence 

The Bayesian evidence is the average likelihood over the entire 
prior parameter space of the model: 



£(e c ,eB)P(ec)P(eB)d Jv ecd M e B , (7) 



where and M are the number of cosmological and Bianchi pa- 
rameters respectively. Those models having large areas of prior pa- 
rameter space with high likelihoods will produce high evidence 
values and vice versa. This effectively penalises models with 
excessively large parameter spaces, thus naturally incorporating 
Ockam's razor. 

The results generated in this paper hav e employed the new 
method of nested sampling dSkillind [2004) as implemented in 
a forthcoming publication (Shaw et al. in preperation). This 
method is capable of much higher accuracy t han previous meth- 
ods such as thermo dynamic integration (see e.g. lBeltran et al.l2005l 
iBridges et all2005h . due to a computationally more effecient map- 
ping of the integral in Eqn. [7] to a single dimension by a suitable 
re-parameterisation in terms of the prior mass X. This mass can be 
divided into elements dX — ir(®)d N ® which can be combined 
in any order to give say 



X(X) 



n{&)d N ®, 



(8) 



'£(e)>A 

the prior mass covering all likelihoods above the iso-likelihood 
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curve C = A. We also require the function C(X) to be a singular 
decreasing function (which is trivially satisfied for most posteriors) 
so that using sampled points we can estimate the evidence via the 
integral: 

Z= f C(X)dX. (9) 
Jo 

Via this method we can obtain evidences with an accuracy 10 times 
higher than previous methods for the same number of likelihood 
evaluations. 



3 BIANCHI VII H SIMULATIONS 

Ja ffe implemented solu tions of the geodesic equations developed 
bv lBarrowetal] Jl985h t o compute Bianch i induced temperature 
maps in type Vllh cases. IJaffe et al.l <2005bl) extended these to in- 
clude cosmologies with dark energy allowing an almost flat geom- 
etry, in keeping with the so called concordance cosmology. One of 
us (ANL, in preparation) independently p erformed the same exten- 
sion and found results agreeing well with lJaffe et ail J2005bl) . 

The models are parameterised by six independent quantities; 
the total energy density Qtot = f2m + ^a, the current vorticity cu, h 
which determines the 'tightness' of the spiral and the Euler angles 
a and /3 determining the position of the centre of the spiral on the 
sky with 7 allowing rotation of the pattern about this point 0. In 
addition one must specify the direction of rotation or 'handedness' 
of the spiral. The ratio of the current shear to the Hubble parameter 
(a/H)o (differential expansion) is determined by the combination 
of these parameters: 

/ a \ _ &h£kot ( ui \ 

^h)o ~ V2(i + hy/z(i + 9hy/2(i~n tot ) \h) ' 

Jaffe found a best fit bianchi template T by a \ 2 analysis of the 
observed CMB sky D for a selection of WMAP data variants by 
minimising: 

X = (D-T) T M _1 (D-T) (10) 

under the assumption of an embedded cosmology M given by 
the b est fit 1-year WMAP angular power spectrum dSpergel et al.l 
120031) . There is no contribution from noise in this expression as 
the dominant uncertainty on large scales is from cosmic variance. 
Their result is shown in Fig. [Trl computed in our formalism with 
parameters Q tot = 0.5, cj = 10.38 x 10" 10 , h = 0.15125 and 
position (42°, 28°, 310°) in Euler angles. The primary feature is 
a central core positioned on the so called 'cold spot', assumed to 
be the source of a large area of non-gaussianity. The corrected map 
(see Fig. |2j, produced after removal of this Bianchi 'contamina- 
tion', provides even by cursory visual inspection a more statisti- 
cally homogeneous sky. 

In Section|4]we will perform our analysis on simulated data in 
order to test our parameter recovery and model selection algorithm 
with a Bianchi component of varying amplitude superimposed on a 
background concordance cosmology. In Section|5]we perform this 



2 We adopt the active zyz Euler convention corresponding to the rotation 
of a physical body in a fixed coordinate system about the z, y and 2 axes by 
7, /3 and a respectively. 

3 See http://www.mrao. cam. ac.uk~jdm57 for an animation of this tem- 
plate on the sphere. 
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Figure 1. Best fit bianchi template as found bv |jaffeetalJl2006l) computed 
from parameters fltot = 0.5, uj = 10.38 X 10 -10 , h = 0.15125 and 
position (42°, 28°, 310°) in Euler angles. 




Figure 2. Raw LILC map (left) with the Bianchi corrected map (right) 



analysis on real data to effectively extend Jaffe's study to fit si- 
multaneously for a selection of right and left handed Bianchi mod- 
els and a background cosmological component. This method has 
two main advantages over Jaffe's original template analysis: i) con- 
firmation of the robustness of the detection within a varying cos- 
mology -thus accounting for possible degenerate parameters and 
ii) estimation of the ratio of Bayesian evidences with and without 
a Bianchi component -providing an effective model selection to 
favour or disfavour the inclusion of a such an addition to the stan- 
dard model. 



4 SIMULATED DATA 

A set of simulated CMB sky maps were produced with the 
HEALPixQ software routines. These contained a cosmic variance 
limited concordance cosmological component and a Bianchi Vllh 
component with fixed handedness, shear, matter and dark energy 
density and geometry but varying the vorticity u> to study the sen- 
sitivity of the method to the total amplitude of the Bianchi signal 
(of which w is a good tracer). Evidence estimates were produced 
for two parameterisations: Bianchi + cosmology (model A); and 
cosmology only (model B). 

In order to provide the best possible assessment of our method 
the Bianchi simulations were created with Jaffe's best fit parameter 
set of: h = 0.15, Q. m = 0.5, Ov = 0, a = 42°, (3 = 28° and 7 = 
310° (see Fig. |3](a)) combined with a best fit WMAP cosmology 
(see Fig. [3] (b)) with eight equally spaced vorticities incremented 
between 5 and 19 x 10 -10 (see Fig.[3](c)-(j)). Ultimately we aim 
to ascertain via these simulations the lowest amplitude of Bianchi 
signal we can succesfully extract via our method -where u — 10 x 

4 http://healpix.jpl.nasa.gov 
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(a) Bianchi 



(b) Concordance cosmology 





(c) oj = 5 X 10" 



(d) oj = 7 X 10" 




(g) u = 13 x 10 



(h) cj = 15x 10" 
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Figure 4. Comparison of the relative amplitudes of the angular power spec- 
tra (C'i) of four simulated Bianchi skies of varying vorticity and the Jaffe 
template. 
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Figure 3. Simulated data created with a WMAP concordance cosmology 
(a) and a Bianchi component (enhanced by a factor of 10 for clarity) (b) of 
the type found by Jaffe with an increasing vorticity (u>). 



10" 



represents the level of sensitivity required to detect a Jaffe 



type component in real data. 

Marginalised parameter constraints (see Fig. [5} illustrate we 
are able to successfully extract the input Bianchi parameters for 
amplitudes at and above roughly u> = 11 X 10 -10 , which suggests 
our method should be able to determine parameters of any Bianchi 
component similar to the Jaffe template, should one exist. As for 
the question of whether such an inclusion is actually warranted by 
th e data, we consi der the evidence. A useful guide has been given 
by I Jeffreys! ( fl96lh to rank relative evidence differences in order of 
significance: AlniS < 1 inconclusive, 1 < AlnE < 2.5 signif- 
icant, 2.5 < AlnE < 5 strong and AlnE 1 > 5 decisive. While 
there is marginal evidence in these simulations for a component at 
u/ = 11 X 10 -10 (see Table [TJ a significant model selection can 
only be made at lo = 13 X 10 -10 , also roughly at the point where 
the Bianchi template becomes visibly apparent in the sky maps in 

Fig.m 




Figure 5. Extracted cosmological and Bianchi parameters from a simulated 
map with ui = 11 X 10~ 10 (described in text). 



5 APPLICATION TO REAL DATA 

The analysis was performed on both one and three year ILC maps, 
probing structure up to multipoles of 64, which does not limit any 
study of Bianchi structure, which contains few features beyond 
I — 20 (see Fig. |4j but does preclude a full examination of the 
background cosmology. This was an unfortunate necessity of using 
the ILC maps owing to their complex noise properties. 

As we are adopting a Bayesian approach in this analysis our 
results will be dependent on the prior parameter ranges chosen. Our 
most conservative, and thus widest uniform prior ranges are listed 
in Table|2] From a Bayesian point of view these are broad enough 
to represent the situation in which there is little prior knowledge 
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Table 1. Differences of In evidence with (model A) and without (model B) 
a Bianchi component with vorticity incremented between 5 — 19 X 10 -10 . 



t<;[xlCr 10 ] A B 



5 


0.0 ±0.2 


-1.6 ±0.3 


7 


0.0 ±0.2 


-1.3 ±0.3 


9 


0.0 ±0.2 


-1.1 ±0.3 


11 


0.0 ±0.2 


-0.9 ± 0.3 


13 


0.0 ±0.2 


+1.4 ±0.3 


15 


0.0 ±0.2 


+3.6 ± 0.3 


17 


0.0 ±0.2 


+6.2 ± 0.3 


19 


0.0 ±0.2 


+8.4 ± 0.3 



Table 3. Cosmological and Bianchi parameterisations for each of the pa- 
rameter subsets studied. 



Model 


Cosmology 


Bianchi 


A 


A s 




B 


A 3 


ng,ng t ,h,u, a, p, 7 ,L/R 


C 


A s 


Qg, Ug t , h, w, a-[0.4,l], /H0.1, 0.7], 7, L/R 


D 


As 


SlB,Jif ot ,fc, W ,7,UR 


E 


As 


, h, w, a, P, 7, L/R 


F 


A s 


f2^, h, ui, a-[0.4,l], /HOT, 0.7], 7, L/R 


G 


As 


f2^, h, u), 7, L/R 



Table 2. Summary of Bianchi Vllh component parameters and priors. 



Full Bianchi 

Qf ot = [0.01, 0.99] 
= [0.01,0.99] 
h = [0.01,1] 
£j = [0,20] x lO - " 
a = [0, 27r]rads 
fi = [0,7r]rads 
7 = [0, 2-7r]rads 
Chirality = L/R 



of the Bianchi models. In particular the h parameter, which de- 
creases the spiral 'tightness' almost exponentially between 0.01 
and 0.2 but little thereafter was given a prior range from 0.01 to 
1 -encompassing the majority of possible features. The vorticity ui 
increases the Bianchi amplitude roughly linearly so a prior range 
centred on the amplitude of the Jaffe template was adopted provid- 
ing a realistic upper limit of 20 x 10~ 10 at which point the Bianchi 
signal would dominate the CMB sky and a lower limit of uj — 
corresponding to no Bianchi signal. Bianchi Vllh models are only 
defined in open cosmologies where Q. to t < 10 so the prior chosen 
on the combination of matter and dark energy energy density was 
restricted to lie within 0.01 and 0.99 -thus introducing an implicit 
prior on of ±0.8 (as + = Qfot)- The Euler angles po- 
sition the centre of the spiral, via a and (3 varied over 360° and 
180° respectively while the orientation of the pattern is given by 
7 rotating over 360° . In the subsequent sections we will analyse a 
range ofBianchi models (Table [3} defined by subsets of the above 
prior ranges while simultaneously fitting a scalar perturbation am- 
plitude A s of an assumed background concordance cosmology. In 
Section[5]2]we will compare, via the evidence, models B - G with 
the concordance cosmology model A as to their ability to describe, 
adequately the observed CMB sky. 

5.1 Parameter Constraints 

The findings of Ijaffe et alj d2005ij) . Ijaffe et all d2005bl) and 
Ijaffe et al] d2006j) all point to a left-handed Bianchi Vllh com- 
ponent with a total energy density Qf ot = 0.5, this value be- 
ing required to focus the centre of the spiral, thus alleviating, at 
least partially, the observable non-Gaussianity dJaffe et ai] [2005a; 

5 Similarly Bianchi VIIo and IX are only defined in flat and closed cos- 
mologies respectively. 
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Figure 6. Comparison of the f2j^ — Bianchi degeneracy (shaded with 
1 and 2a contours) with the familiar CMB geometric degeneracy from 
WMAP first (blue) and third year + polarisation (green) data (with 1 and 
2cr contours). 

iMcEwen et al]2005l ; lMcEwen et al]2006|) if centered in the south- 
ern hemisphere. However both 1-year and 3-year WMAP data 
would suggest a value of Qtot much closer to unity, perhaps larger, 
motivating the search for a Bianchi model within a more viable cos- 
mology. The extension of the Bianchi parameter space to include a 
dark energy densi ty fi^ should thus b e explored thoroughly to alle- 
viate this tension. Ijaffe et all (2005b) found that a degeneracy is in- 
troduced with fi^ in the Q.^ — Q.^ plane, similar to the 'geometric' 
degeneracy which exists in the CMB. Along this region, Bianchi 
skies appear identical, only being distinguished by their relative 
amplitudes. The maximum ln-likelihood of a Bianchi component 
at = 0.5, f^A = is -2716.31. At locations that are broadly 
consistent with concordance, such as i!^ = 0.3, = 0.69 or 

= 0.2, ttf = 0.79 this falls to -2717.81 and -2719.11 respec- 
tively. A more thorough exploration of the space by MCMC sam- 
pling (Fig.[6]) confirms that despite the increased degree of freedom 
that Sl^ introduces, at the lcr level no overlap with either 1 or 3- 
year WMAP likelihood contours is found. Based on these results 
we can effectively rule out Bianchi Vllh models as the physical 
origin of this effect. Of course the corrections that a Bianchi com- 
ponent makes to the WMAP data are no less interesting, whether 
we can model them or not. This is especially true given that the ar- 
eas of non-Gaus sianity have been mo stly preserved into the 3-year 
WMAP release dMcEwen et alj2006t) . 

Before continuing it is worth discussing in some detail the ef- 
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feet of each parameter on the structure and scale of the Bianchi 
templates. The morphology is most sensitive to and h, which 
determine the 'focus' and 'tightness' (also to some degree the am- 
plitude) of the spiral, respectively. The introduction of dark energy 
in the form of a cosmological constant energy density creates a 
correlation between the matter density as the repulsive nature of the 
former de-focuses the spiral. This degeneracy is almost exact, being 
broken only by the overall amplitude. Unfortunately, many other 
parameters can influence the amplitude of the CMB anisotropics on 
this scale such as the amplitude of primordial scalar (and tensor) 
density perturbations, the optical depth to reionisation, the spec- 
tral 'tilt' of the primordial perturbations and of course the Bianchi 
vorticity uj. This profusion of possible (and mostly degenerate) ex- 
planations for the CMB amplitude on large scales make constraints 
difficult to obtain. Given these difficulties we feel it is justified to fix 
the cosmological pa rameters at their WMAP 3-year best fit values 
(Spergel et al. 2003) varying only the scalar amplitude A s , More- 
over, careful choice needs to be made over the priors placed on the 
Euler angles particularly as large regions of the ILC maps contain 
residual Galactic noise and other emission. In this light we believe 
a reasonable restriction on the position (a and (3) of the centre of 
the pattern is well motivated, while leaving the orientation angle 
(7) free to rotate the spiral. Parameter constraints with and without 
Qf (models C and F) are shown in Figs [7] and [8] for a left-handed 
Bianchi component. The angular ranges used, a — [0.4, 1.0] and 
(3 = [0.1, 0.7] do implicitly place a prior on the centre of the fea- 
ture being at least partially associated with the 'cold spot' in the 
southern hemisphere located at roughly a — 0.7, f) = 0.4. Right- 
handed models examined were entirely unconstrained by the data 
(see Fig. 9) with a typical best fit In likelihood at l east 2 units below 
left-handed models, substantiating the claim in Ijaffeet all (2006) 
that the only viable Bianchi fit is left-handed. The degenerate ef- 
fect of is clearly evident in the laxer constraints, and since we 
are now only assessing the feature as a non-physical effect there is 
no benefit in retaining the parameter -it provides no additional ef- 
fect on the template structure. The above constraints pro vide robust 
confirm ation of a template very similar to that found by jjaffe et all 
l2005al) . 



5.2 Model Selection 

The conclusion that the Bianchi Vllh models may only be treated 
as templates does not preclude a consideration of the necessity of 
such a template to explain the data. Indeed, removal of areas of 
non-Gaussianity and correction of large-scale power would explain 
serious anomalies in the current data, but at the cost of removal of 
universal isotropy. The Bayesian evidence provides an ideal way to 
assess quantitatively whether such action is necessary. As with the 
previous section we shall assess a template both with and without a 
dark energy density and with a range of priors on the two positional 
Euler angles via models B - G. All model comparisons are depicted 
via In evidence differences; hiA — lne-G so that a positive value 
illustrates a preference for the inclusion of the Bianchi component 
and vice versa. 



In agreement with the parameter constraints, little evidence 
exists for any of the right-handed models, typically at least 1 In 
evidence unit lower than for the left-handed models. We found the 
results were relatively insensitive to the choice of priors on the en- 
ergy densities, h and ui, but heavily dependant on the Euler angles. 




0.2 0.4 B 0.6 0.8 



tot 

Figure 7. Extracted Bianchi Vllh parameters including (left handed 
model C) from WMAP 3-year data. 




Figure 8. Extracted Bianchi VII^ parameters with a prior of Sl^ = (left 
handed model F) from WMAP 3-year data. 



The only preferred model (a marginally significant In evidence dif- 
ference > 1) was found by restricting the pattern to lie centred at 
the 'cold-spot' (models D and G), parameter constraints are shown 
in Fig. [TO] It must be noted that the Bayesian evidence is a prior 
dependent quantity, so that in general a smaller prior will tend to 
lift the evidence and vice versa, thus an increased detection using 
these restricted priors is understandable. However, by still allow- 
ing the pattern to rotate freely, via 7 we still gain all the degrees 
of freedom necessary to allow the data to decide the best morphol- 
ogy. Indeed all we are saying with this prior is that we believe the 
'cold-spot' to be driving any detection. The effect of increasing the 
parameter space with the addition of is clealy offset by the large 
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Table 4. In evidence differences using one and three year all sky maps ILC maps for L and R handed models including a dark energy density. 



Data \ Model 


B(L) 


B(R) 


C(L) 


C(R) 


D(L) 


D(R) 


3-year ILC 
1 -year ILC 


-0.9± 0.2 
-0.9± 0.2 


-1.4 ±0.2 
-1.5 ±0.2 


-0.6 ± 0.2 
-0.5 ± 0.2 


-1.2 ±0.2 
-1.1 ±0.2 


+1.2 ±0.2 
+1.2 ±0.2 


-1.0 ±0.2 
-1.0 ±0.2 


Data \ Model 


E(L) 


E(R) 


F(L) 


F(R) 


G(L) 


G(R) 


3-year ILC 
1-year ILC 


-0.9± 0.2 
-0.8± 0.2 


-1.3 ±0.2 
-1.5 ±0.2 


-0.6 ± 0.2 
-0.7 ± 0.2 


-1.1 ±0.2 
-1.1 ±0.2 


+1.1 ±0.2 
+0.9 ± 0.2 


-1.0 ±0.2 
-1.0 ±0.2 



6 CONCLUSIONS 

We have investigated the possibility of signatures of universal shear 
and vorticity in 1 and 3-year WMAP data. We implemented Bianchi 
Vllh simulations in a fully Bayesian MCMC analysis to extract the 
best fitting parameters 0. of the model and to determine whether its 
inclusion is warranted by the data. We only make a marginally sig- 
nificant detection for a left-handed Bianchi template with parame- 
ters; Qtot = 0.43 ± 0.04, h = 0.32±g;?§, w = 9.7±i;| x 10" 10 
with orientation 7 = 337° I21) if we restrict its position on the sky 
to lie centred on the 'cold-spot' at (a = 42°, (3 = 32°) -all other 
models are essentially disfavoured by the data in both 1 and 3 year 
maps. Since the open universal geometry implied by these tem- 
plates effectively rules out the Bianchi Vllh models as the origin 
of this structure we can only treat these models as template fittings, 
morphologically similar to some as yet unknown physical cause. 
Nevertheless our analysis demonstrated that there are interesting 
effects on large scales yet to be explained in the WMAP data. 




1 2 3 4 6 8 0- 2 0.4 #.6 0.8 



P[rad] y[rad] " m 

Figure 9. Extracted Bianchi Vllh parameters with a prior of f!^ = (right 
handed model B) from WMAP 3-year data. 




u.o u.o " 1n 

h o[x10- 10 ] 




Figure 10. Extracted Bianchi Vllh parameters with a prior of = and 
centered on the 'cold spot' (left handed model G) from WMAP 3-year data. 
The only model considered which yields a positive evidence. 



volumes of high likelihood that exist along the degeneracy, leaving 
the evidences almost unchanged. The conclusions are essentially 
independent on going from 1 to 3 year data, presumably as both 
maps are effectively cosmic variance limited up to at least I < 10 
where the majority of Bianchi structure lies. 
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APPENDIX A: BIANCHI VII H TEMPLATES 



We illustrate here our best fitting Bianchi VHh templates found 
from 3 year WMAP data. 



